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ABSTRACT 

Numerous X-ray bursts show strong oscillations in their flux at sev- 
eral hundred Hz as revealed by RXTE. Analyzing one such oscillation 
from the X-ray binary Aql X-1, I find that low energy photons (3.5-5.7 
keV) lag high energy photons (> 5.7 keV) by approximately 1 radian. 
The oscillations are thought to be produced by hot spots on the spinning 
neutron star. The lags can then be explained by a Doppler shifting of 
emission from the hot spots; higher energy photons being emitted earlier 
in the spin phase as the spot approaches the observer. A quantitative 
test of this simple model shows a remarkable agreement with the data. 
Similar low energy lags have been measured in kilohertz quasi-periodic 
oscillations and in the accreting millisecond pulsar SAX J1808. 4-3658. 
A Doppler delay mechanism may be at work there as well. 



Subject headings: accretion, accretion disks — black holes - stars: neu- 
tron — X-rays: stars 
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1. Introduction 

The Rossi X-Ray Timing Explorer (RXTE) 
has uncovered strong oscillations of X-ray flux 
during X-ray bursts in several low mass X- 
ray binaries (Strohmayer et al. 1996). Cur- 
rent interpretation favors a rotation mecha- 
nism for the burst oscillations: asymmetric 
nuclear burning leaves a 'hot spot' which ro- 
tates with the neutron star and produces a 
strong modulation (Strohmayer et al. 1998). 
The frequency of the burst oscillation is then 
the spin frequency of the neutron star, or 
twice the spin frequency for two spots (Miller 
1999). Oscillations have been discovered in 
X-ray bursts from the following systems: 4U 
1728-34 (363 Hz, Strohmayer et al. 1996; 
Strohmayer, Zhang & Swank 1997), KS 1731- 
260 (524 Hz, Smith, Morgan & Bradt 1996), 
a source near the galactic center (589 Hz, 
Strohmayer, Jahoda & Giles 1997), Aql X- 
1 (549 Hz, Zhang et al. 1998), 4U 1636- 
536 (581/290 Hz, Strohmayer et al. 1998, 
Miller 1999), and 4U 1702-429 (330 Hz, Mark- 
wardt, Strohmayer & Swank 1999). The ob- 
served frequencies are close to the 401 Hz spin 
frequency of the accreting millisecond pul- 
sar SAX J1808.4-3658 (Wijnands & van der 
Klis 1998), further strengthening the identifl- 
cation of these frequencies with the neutron 
star spin. 

The detailed energy dependence of these 
burst oscillations is one avenue that remains 
to be explored. Here I show that the low en- 
ergy photons in a burst oscillation from Aql 
X-1 lag the high energy photons by roughly 
15% of the oscillation period. Lags of the 
same sign and similar magnitudes have also 
been detected in other fast signals from low 
mass X-ray binaries: the kilohertz quasi-periodic 
oscillations, QPOs (Vaughan et al. 1997, 
1998, Kaaret et al. 1999) and the SAX 



J1808.4-3658 pulsed emission (Cui, Morgan 
& Titarchuk 1998). 

A simple mechanism of Doppler shifted 
emission may explain these lags. Strong Doppler 
effects are expected to be important since the 
fast spin rates imply high speeds [fi — vjc 

0. 1.. As a hot spot on the spinning neu- 
tron star approaches the observer (at early 
phases) the emission is Doppler boosted and 
blue shifted, as it recedes (at later phases) 
the emission is deboosted and red shifted. At 
early phases the spectra are also attenuated 
due to the smaller projected area. The re- 
sult is that low energy photons arc preferen- 
tially emitted after the high energy photons. 
A quantitative test of this Doppler delay sce- 
nario matches the observed low energy lags 
in Aql X-1 well. The possibility of Doppler 
effects and the fact that they may manifest 
in pulse phase spectroscopy has been noted 
before by Strohmayer et al. (1998). 

In the next section I present the measure- 
ment of the lag in the X-ray burst from Aql X- 

1. In Section 3 1 describe a simple model for 
the relativistic effects and compare the pre- 
dicted delays to those observed. Section 4 
discusses these results in a broader context. 

2. MeEisurements 

For this analysis I consider the X-ray burst 
from Aql X-1 starting 1 1997 March 1 23:27:40 
UTC (see Zhang et al. 1998 for a report of 
this burst). I use data from the RXTE Pro- 
portional Counter Array (PCA) in an 'event' 
mode with high time resolution (122//sec) and 
high energy resolution (64 channels). A sec- 
tion of the lightcurve is shown in Figure 1 
(top). There are gaps in the event mode 
data since the required telemetry rate is high. 
Within the 4 second time window shown in 
Figure 1 (top), the power density spectrum 
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for all the channels shows a strong oscillation 
at 549.7 Hz (Figure 1, bottom). In the follow- 
ing 1 calculate Fourier transforms within this 
time window. 

Phase delays in a signal between two en- 
ergy bands are quantified by means of cross 
spectral analysis (van der Klis et al. 1987; for 
more information see Vaughan et al. 1994; 
Nowak et al. 1999). The cross spectrum is 
defined as C(j) = X^(j)X2(j), where X are 
the measured complex Fourier coefficients for 
the two energy bands at a frequency Vj. The 
phase lag between the signals in the two bands 
is given by the argument of C (its position an- 
gle in the complex plane). The error in the 
phase lag is calculated here from the coher- 
ence function uncorrected for counting statis- 
tics (Nowak et al. 1999). The cross corre- 
lation code used here has been employed to 
calculate phase lags in black hole candidates 
(Ford et al. 1999), SAX J1808.4-3658 and 
kilohertz QPOs and matches the results re- 
ported in the literature. 

Figure 2 shows the resulting phase lags 
from the cross spectra of the 4 seconds of 
data described above. Negative numbers in- 
dicate that the oscillations in the low energy 
band (3.5-5.7 keV) lag those in the higher 
energy bands. The lags are calculated by av- 
eraging the signal in the range 549.6-550.1 
Hz. The delays in each band up to 30 keV 
(where background dominates) are 3(T signifi- 
cant. The delay between 3.5-5.7 keV and the 
entire 5.7-43.6 keV band is 0.93 ± 0.18 rad, 
5(T significant. 

Deadtime effects can in principle affect the 
measured phase lag. The data considered 
here are in the tail of the burst (rate of 9280 

c/s, full energy band) where deadtime is less 
important. One method of correcting for 
deadtime is to subtract a cross vector aver- 



aged over high frequencies where no correla- 
tion is expected (van der Klis et al. 1987). 
Employing this correction does not change 
the values measured here. 

Due to the data gaps it is not possible to 
perform cross correlations on long stretches 
of data earlier in the burst. Cross correla- 
tions on 0.5 sec intervals of data earlier in the 
burst return large errors on phase delays with 
inconclusive results. 

3. Model 

As a simple model for the lags 1 consider 
discrete hot spots on the surface of the ro- 
tating neutron star. The rest frame emission 
of the clump is a blackbody. The observed 
spectrum at frequency v at spin phase is: 

F^{v) = Aocos5[7-^(1 - (3iicose)-^f 

where (3 = v/c, 7 is the Lorentz factor, 
kT = kTQj~^{l—i3^cos9)~^ (with kTo the rest 
frame temperature), is the sine of the an- 
gle between the spin axis and the line of sight, 
and Aq is a normalization. The above formula 
is a relativistic transformation of the black- 
body which shifts kT and modifies the nor- 
mahzation such that F/u^ is conserved (see 
Rybicki & Lightman 1979). The cos5 term 
is an area projection factor, with S the an- 
gle between the normal and line of sight in 
the rest frame {5 ^ tt — 9). The phase angle, 
9, is defined such that phase zero is with the 
spot approaching the observer directly. The 
spots are considered small and isotropically 
emitting in the rest frame. 

1 take kTo=l keV, (3 = 0.1. These are val- 
ues appropriate for the neutron star; a more 
exact value of kTo is in principle possible from 
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the spectral fits but this depends on the frac- 
tion of the surface contributing to the modu- 
lated hot spot emission. A more exact value 
of /5 depends on the neutron star radius. I 
also take /i — 1, i.e. a line of sight through 
the equator. The spin frequency is 275 Hz 
and two antipodal hot spots produce an os- 
cillation at 550 Hz. Such a geometry, where 
the ~550 Hz signal is a harmonic of the spin, 
is suggested by recent results on other burst 
oscillations (Miller 1999). The resulting spec- 
tra are blackbodies whose temperature shifts 
by 10% over the period. Averaged over phase, 
the spectrum is approximately blackbody in 
shape with kT within 1% of the input kTo. 

The spectra as a function of 9, folded 
through the RXTE response matrix, yield 
lightcurvcs of count rates in various energy 
bands. From these lightcurvcs 1 calculate 
the phase lag in the 550 Hz signal with the 
FFT and cross-correlation programs used in 
the measurements above. The results of this 
calculation arc shown with the data in Fig- 
ure 2. There arc no free parameters, only 
the assumptions taken above. The calcu- 
lated lags will decrease if kTo is increased or 
P is decreased. The smaller delay for higher 
kTo happens since the peak of the lightcurvcs 
comes later in phase for higher energy pho- 
tons, corresponding to a smaller delay be- 
tween high and low energy photons. Observ- 
ing at higher inchnations (decreased //) will 
also decrease the lag. The lightcurvcs that 
yield these predicted lags generally have max- 
ima at earlier phases for higher energies and 
are more sharply peaked in shape at higher 
energies. 

This simple model neglects general rela- 
tivistic effects (e.g. Strohmayer 1992, Miller 
& Lamb 1998). Two main factors from GR 
will effect the observed light curves. Gravita- 



tional bending makes the spots observable at 
^ < or ^ > TT, stretching the pulse. Light 
travel time delays, longer for more extreme 
bending, will also shift the pulse. These ef- 
fects depend on the compactness of the star. 
Given the quality of the present data, a more 
detailed treatment including these effects is 
not justified. An overall gravitational red- 
shift also means that kT in the local frame 
is higher, as in X-ray burst spectral models. 

4. Discussion 

The previous sections show that low energy 
photons lag high energy photons in the oscil- 
lation signal of an X-ray burst from Aql X-1. 
The sign and magnitude of the lags are in 
agreement with the simple model considered 
in Section 3 of two hot spots on the neutron 
star producing a Doppler boosted and shifted 
spectra as the star rotates. 

This Doppler delay mechanism for produc- 
ing low energy lags may describe not only the 
lags in the X-ray burst oscillations but also 
the lags in the accreting millisecond pulsar 
SAX J1808.4-3658 (Cui, Morgan & Titarchuk 

1998) and the (lower frequency) kilohertz QPOs 
(Vaughan et al. 1997, 1998; Kaaret et al. 

1999) . Both show a lag of low energy photons 
relative to high energy photons with mag- 
nitudes of roughly ~100 /isec (~ 0.3 rad) 
for SAX J1808.4-3658 and ~30 ^usec (~ 0.2 
rad) for the kilohertz QPOs in similar energy 
bands to those considered here. Some models 
link the frequency of the kilohertz QPOs to 
a Keplerian motion in the disk (Miller, Lamb 
& Psaltis 1998, Stella & Vietri 1999; but see 
Titarchuk, Lapidus & Muslimov 1998). If any 
of the kilohertz QPOs is a result of Keple- 
rian motion, one might expect a soft lag due 
to Doppler delays. Such lags have been ob- 
served in what is likely the lower frequency of 
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the two QPOs. 

Doppler delays are an alternative to pre- 
vious mechanisms invoked to produce lags. 
Comptonization has been one process used to 
explain low energy lags in SAX J1808. 4-3658 
(Cui, Morgan & Titarchuk 1998). Low energy 
lags are produced if high energy photons are 
injected into a relatively cool Comptonizing 
cloud. This is the opposite of the situation 
normally considered: Comptonization by a 
hot cloud in the same region. A hot cloud pro- 
duces a lag of high energy photons, as shown 
quantitatively for fast signals by Lee & Miller 
(1998). Another mechanism suggested for low 
energy delays is an extended, cooling hot spot 
with lower energy photons from the outer re- 
gions (Cui, Morgan & Titarchuk 1998). 

More measurements of phase lags in X-ray 
burst oscillations are clearly needed, in par- 
ticular in the ~350 Hz oscillations which are 
likely from single spots. Improved statistics 
will also yield a better test of the predicted 
energy dependence of the lags. 
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Fig. 1. — The lightcurve (top) and power spectrum (bottom) for the 1997 March 1 X-ray burst 
from Aql X-1. The power spectrum is from four seconds of data taken in the time window shown 
by the dashed hnes. Data gaps are due to telemetry overloads in this PCA mode. 
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Fig. 2. — Phase delay measurements for the X-ray burst oscillation shown in Figure 1 relative 
to the 3.7-5.7 keV band. A negative value indicates that low energy photons lag high energy 
photons. The solid line is the Doppler delay model for two hot spots (see text). 
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